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Abstract: We consider light trapping in an amplifying medium consisting of cold alkali-metal 
atoms; the atomic gas plays a dual role as a scattering and as a gain medium. We perform Monte- 
Carlo simulations for the combined processes. In some configurations of the inhomogeneous distri¬ 
bution this leads to a point of instability behavior and a signature of random lasing in a cold atomic 
gas. 


Introduction 

Radiation trapping [l[, light localization [2] and cooperative phenomena involving coherent matter-radiation 
couplings are the subjects of various and intensive discussions in the scientific literature. Among areas of current ex¬ 
amination are exploration of super- and sub-radiant configurations M , cooperative effects including the cooperative 
Lamb shift (fUlljj. and the role of collective phenomena in dense and cold atomic gases [mil- 

In this paper, we concerned mainly the intriguing phenomenon of the so-called random laser Il7i - l20| predicted 
originally by V. Letokhov in 1968. As was shown in his seminal paper E3, if the amplifying active medium is 

distributed throughout another elastically scattering medium, and if the sample size r o fulfills the threshold condition 

1 /2 

7*0 > 7T (Itrlg /3) ' , where Z tr is the transport length of light diffusion and l g is the gain length, there is an instability 
point in the amplification process. Of course, in a realistic scenario the infinite amplification is impossible. If the 
sample volume is large enough such that this threshold is overcome then the amplification initiates a saturating 
feedback of radiation with the active medium and finally the radiation of the entire system can approach a steady 
state regime. The role of randomly distributed elastic scatterers in this process consists of effective trapping of the 
amplified radiation such that the scatterers serve as the cavity in a conventional laser scheme. 

There have been many studies of random lasing, and signatures of the threshold behavior of the random lasing 
phenomenon in condensed matter systems. For instance, a narrow spectral feature was observed in the emission of 
metal ions distributed in powder paint [TO ] . Another example is the radiation of a dye active medium, which can be 
diffusely trapped by TiC >2 microparticles distributed in the dye solution [22t | . These and other examples of laboratory 
manifestations of the phenomenon are reviewed in fl8j . The physical evidence of lasing from the disordered system is 
probably also observable in astrophysics as stimulated emission in certain stellar objects [HI . In laboratory realizations 
it seems of special interest to systems with a quasi one dimensional geometry where the lasing in a disordered medium 
can be accompanied by localization phenomenon f2~i1 . [ 2 I ] . 

However, in the gas phase (e.g. atomic or molecular gases) unambiguous demonstration of all expected features of 
random lasing is still a challenging problem for experimental verification. It has been proposed that this effect could 
be observed in a system of cold alkali-metal atoms prepared in a magneto optical trap [26(. The unique statement 
formulated in that paper is that the radiation can be emitted and trapped by the same medium, which is active and 
scattering simultaneously. Such a system can consist of billions atoms at the temperature ~ [iK with very large 
optical thickness. The authors made their estimates based on a two-level model of optical transition pumped by a 
strong coherent field, i.e., for the Mollow-type system and they predicted an optical depth around several hundreds 
to approach the threshold condition. But as was recognized later, the complicated multilevel structure of an alkali- 
metal atom would be even more helpful in realization of the random lasing conditions. Indeed in the alkali-metal atom 
system the atoms can be repopulated among different hyperfine sublevels and then can be considered as distinguishable 
fractions of the matter. The progress in architectures of cold atomic systems is rapidly developing and the required 
conditions under a variety of circumstances may be attainable in the laboratory in the near future. Recently, elegant 
experimental results showing that the threshold conditions can be overcome in the system of cold atoms have been 
reported in [23, [HI ■ 
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In the present paper we discuss a particular and different scheme showing how both amplification and light trapping 
can be organized in such a system. For this scheme, we show that threshold conditions can be attained, leading to a 
regime of random laser generation. We also briefly discuss how one component of such a scheme, an optically induced 
and quasi one dimensional soft cavity may be generated. 


Amplification of the spontaneous Raman emission via 
radiation trapping mechanism 

In this section we consider the mechanism of Raman-type amplification for radiation trapped in a disordered 
atomic medium. Such a mechanism can be organized in the Z? 2 -line hyperfine manifold of any alkali-metal atom and 
in Fig. [T] we show 85 Rb as a relevant example. The respective excitation geometry is reproduced in Fig. [2] Let the 
rubidium atoms initially populate the upper hyperfine sublevel in the ground state. Then some portion of them can 
be transferred to the lower sublevel with a microwave field. Upon applying a strong coherent field near resonance 
with the electric-dipole forbidden Fo = 2 —» F = 4 transition one can initiate the off resonant Raman process via 
F = 3, 2 upper states, thus transferring atoms back to the Fo = 3 sublevel. The radiation so generated will have a 
frequency nearly degenerate with the practically closed Fo = 3 —> F = 4 transition. The crucial peculiarity of such 
an excitation scheme is that the light emitted in the spontaneous Raman process can be strongly trapped by elastic 
scattering on the closed Fo = 3 —» F = 4 transition if the optical depth of the sample is large. Then the portion of 
light, originally created in the spontaneous Raman process and then propagating via multiple scattering through the 
disordered atomic ensemble can stimulate additional Raman emission. If we imbalanced the trapping conditions in 
such way that the amplification of the transported light overcome its losses, then we could arrive at the situation of 
a random lasing mechanism as discussed in [291 ]. We briefly summarize that result below. 

The above conditions justify that only a small portion of the atoms are repopulated from the upper hyperfine 
sublevel of the ground state and the stimulated amplification is expected to take place at a longer path length than 
the diffusion transport length. Then the amplification process can be subsequently considered in the perturbation 
theory and conveniently described in the diagram framework beyond any simplifications associated with light transport 
or diffusion equation. The electric field correlation function can be given by the sum 

/ OO . . 

exp(-Mwr) T)E { +\r u ty) dt , 

-OO ' ' 

d[\] m = M + £>$ M + D§ H + ... (0.1) 

where are the Heisenberg operators of the positive/negative frequency components of the electric field. For 

the sake of brevity we incorporated all the notation details into the numerical arguments such that 1 = /n, n,... etc. 
and indicate by the numbers the different arguments. At each n-th scattering step the dynamics of the field obeys the 
following Bethe-Salpeter-type transformation for the correlation function of light accumulated in the previous n — 1 
scattering sequence: 

^ii' ( w ) ^ 'y ' y ] Pm 3 ,m3 X 

1712=1712 m3,771,3 

GF 0 =2,3 eF 0 =3 

v[y ] (w) *a 23 (w) * {u)*a\, 2 ,{uj)* V$, (w), (0.2) 

where asterisks denote here the integral over spatial (if necessary) and sum over polarization and internal atomic 
variables. The zero order is associated with the light originally emitted in the spontaneous Raman process. Here 
T>[tp (w) and T>[^ (w) are the retarded and the advanced-type photon propagators, which express each other by 
Hermitian conjugation such that T>^ = D^)*. The internal part of the transformation is the matrix product of 
the correlation function for light incident on a random atom with the scattering tensor a 23 (w) and its Hermitian 
counterpart u\, 2 ,(u>) weighted with the density matrix of atoms populating the state Fo = 3 (Fig. [T]). The 

sign of proportionality in the transformation (10.21) indicates any required factors associated with relevant scaling and 
geometry. 

The light propagation dynamics mostly depends on the following kinetic parameters responsible for the light trans¬ 
port. The imaginary part of the sample susceptibility determines the extinction length for the plane wave entering 
the sample at the e -1 level of losses. The inverse extinction length is given by 

CV) = n 0 cr ex (u}) = 4nk Imy(w), 


( 0 . 3 ) 
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where no is a typical local atomic density and for the sake of simplicity we omitted in this estimate its spatial 
dependence. In an anisotropic sample this quantity can be defined only for a plane wave propagating along a specific 
direction associated with the main reference frame, such that x( w ) i n Eq. (10.31) can be any of the major components 
of the susceptibility tensor. Otherwise the definition has to depend on the propagation direction and track the 
polarization properties of light. It is important to recognize that in the case of amplification the extinction length 
accumulates not only losses but also the gain associated with the stimulated Raman scattering. This parameter can 
be even negative in the case of population inversion and then indicates the distance of e +1 amplification. 

Another kinetic parameter responsible for the scattering process, and in particular for light trapping, is the scattering 
length, which is given by 

CcV) 

O'sc(w) 

where we introduced the scattering cross section a sc (co ) with the scattering tensor, which for the scattering on a 
dressed atom is specified in (29j |. Here w, e and cj', e' are the frequencies and polarization vectors of the input and 
output photons respectively and the integral is over the full scattering angle. Similarly to the extinction length, this 
quantity critically depends on the polarization direction of the incident photon. 

Light diffusion mostly depends on the elastic contribution when w = u/ and both m and m! belong the same Fq = 3 
level. With keeping only elastic contribution we can define the characteristic length associated with losses, 

CM = -CM = CM - CM, (0-5) 

which evaluates the averaged distance that a photon travels before being lost by an inelastic scattering event. The 
first equality indicates that in an amplifying medium it is possible to have this quantity negative in sign and then 
redefine it as the gain length l s (u>). In such a gain medium the radiation trapping and amplification mechanisms can 
overcome the instability point, entering the regime of random Raman laser generation. 

In Fig. 0we show one example of the Monte-Carlo simulations from j29j for the spectral variation of the intensity 
of the light originally emitted by the Raman process and further scattered by the atomic ensemble under the trapping 
conditions. The observation channel is chosen at the angle 9 = 45° to the polarization direction of the driving fields. 
In the presented data the frequency of the optical driving mode was scanned near the forbidden F 0 = 2 —>• F = 4 
transition, for particular Rabi frequencies of the control optical and microwave modes, and the optical thickness of the 
atomic sample on the Fq = 3 —> F = 4 transition was varied from lower to higher values given by b o = 1, 5,10,15, 20. 
For the sake of simplicity we omit in our discussion some specific details concerning the coordination of the microwave 
and light induced energy shifts and the problem with spectral distribution of the emitted light [29| . 

The graphs of Fig.[3]show that in the broad spectral domain the emitted radiation undergoes a sequence of scattering 
on rather complicated Autler-Townes resonance structure created by the driving fields. The susceptibility spectrum, 
modified by both the driving optical and microwave modes, was calculated in [29|. It is crucial that the X^-type 
nonlinearity (with respect to the optical excitation) manifests itself in creation of additional quasi-energy resonances 
located in the atomic spectrum near the frequency of the optical mode as indicated by the dashed bar in Fig. |T| The 
respective quasi-energy level always interacts with the modes emitted in the Raman process and traps them even 
if the optical mode is tuned far off-resonant from the F = 4 reference level, i. e., in normally transparent spectral 
domain. But such trapping mechanism is not closed because the interaction via the quasi-energy level opens inelastic 
inverse anti-Stokes scattering channels as well. The anti-Stokes scattering process redistributes the Raman emission 
into those spectral modes, which further escape the sample without any absorption, and in the context of the random 
lasing effect this process should be treated as an additional loss mechanism. When the optical mode is scanned 
in the vicinity of the upper level F = 4 the situation noticeably changes and the trapping effect, also initiated by 
elastic scattering on the closed F$ = 3 —>• F = 4 transition, leads to additional amplification, which is visualized as 
a bump-shaped enhancement of the spectra. However as follows from the presented data the significant part of the 
light still escapes the sample via inelastic inverse anti-Stokes scattering channel. Although this type of scattering 
essentially reduces the light intensity transported via the elastic channel, when the optical mode scans the vicinity 
of the Fo = 2 —>• F = 4 transition it nevertheless stimulates Raman amplification of the trapped light, which can be 
seen in Fig. [3] in both the elastic and inelastic channels. 
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Raman emission and light trapping in spatially inhomogeneous conditions 


As was commented on in (2f| the problem with the extra losses, related to the control mode nonlinearity, can be 
solved in a spatially inhomogeneous configuration. The crucial requirement for that is to separate the spatial location 
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of the amplification and trapping areas. In the amplification area the atoms have to populate only the pumping 
level, i. e., F 0 = 2 in the considered case. That eliminates any inverse scattering processes. If this active area was 
surrounded by the atoms trapping the emitted radiation that would induce a certain soft cavity feedback and can 
create instability in the system. 

In Fig. [4] we show one possible experimental architecture, which implies preparation of a spatially inhomogeneous 
energy structure and population distribution of the hyperfine sublevels in the atomic ensemble. This can be performed 
via controllable light shift of only one hyperfine energy level for the atoms located in a spatially selected volume of 
cylindrical symmetry inside the atomic cloud. If we assumed the level Fo = 2 as light shifted and organized the 
population inversion for all the atoms in the selected volume onto this level (for example, with 7r-type microwave 
pulse) then these atoms could form an active medium for the photon emission. Indeed, as follows from the transition 
scheme, shown in Fig. [TJ and explaining diagrams in Fig. [2 the control mode would create the photon emission on 
F = 4 —>• Fo = 3 only for the atoms inside the selected volume and the interaction would be off-resonant with the 
control field for the other atoms of the ensemble. Then the atoms outside of the active volume would trap the emitted 
light and play the role of a soft cavity redirecting the light in a quasi-one-dimensional propagation channel associated 
with that volume. That could lead to instability if each spontaneously emitted photon would have a diffusion path 
long enough for stimulation of extra photon emission while it propagates through the channel. 

In Fig. [5] we show an example of our Monte-Carlo simulation performed for such a system under conditions demon¬ 
strating the instability behavior. These graphs track how the intensity of a single external light source located in a 
central point inside such a soft cavity is expanded in different orders of light scattering and depends on the amplitude 
of the pumping initiated by the control mode. At low pump intensity we observe an amplification of the outgoing light 
intensity but converging sum over a limited number of the scattering orders. At high intensity of the control mode the 
process becomes diverging and unstable. This effect depends on the optical density of the surrounding atoms &o and 
is more evident as far as the optical depth is higher. We can associate this instability with transforming the emission 
process towards the regime of random laser generation. 


Creation of a soft cavity 

Preparation of a controlled realization of a soft cavity for this and other potential applications is a challenging 
experimental enterprise. But some basic elements of the soft cavity configuration have been already demonstrated in 
the laboratory |30|, |3l| • One important constraint is that the density of the atomic sample be relatively high, so that 
radiation following optical excitation within the soft cavity has a tendency to be confined within the cavity. This 
suggests that one route towards experimental implementation would be to use an optical dipole trap to achieve the 
necessary high density. A second constraint would be that the soft cavity should be transversely small, on the order 
of a few microns or less, in order to limit the number of transverse optical quasimodes within the approximately 
cylindrical channel geometry. 

We have realized such a quasi-one-dimensional configuration by optical manipulation of a micron-scale optical 
dipole trap (FORT) containing 8 'Rb atoms. The dipole trap is loaded from a magneto optical trap using procedures 
described elsewhere [30j, [H|. The atomic samples so produced have a peak density ~ 2 x 10 13 atoms/cm 3 and a 
characteristic temperature around 100 [iK. A soft cavity is generated by using a combination of optical excitation 
on hyperfine components of the D 2 and Di lines of 87 Rb. Two optical beams are used; the beams are combined in a 
single mode optical fiber, and focused tightly and cross wise to the long axis of the optical dipole trap. One of the 
beams is a control beam, with a peak power on the order of 100 mW, and is tuned in the vicinity (1-100 GHz) of the 
Di transition around 795 nm. We term this beam the light shift laser. The second beam is tuned close to the nearly 
closed Fo = 2 —> F = 3 hyperfine transition associated with the D 2 line at a wavelength around 780 nm. The intensity 
of this probe beam is very low, well lower than the on-resonance saturation parameter for the probed transition. The 
transmission of this beam through the atomic sample serves as a probe of the influence of the light shift laser. In order 
to suppress detection of the light shift laser in the forward direction, several optical interference filters are used; these 
transmit a convenient fraction of the probe beam, while reducing the light shift laser beam intensity to undetectable 
levels. We point out that it is essential to also filter the output of the lasers generating the probe and light shift 
laser beams so that broadband fluorescence from the diode outputs do not influence the measurements. Additional 
interference filters and colored glass filters are used for this purpose. 

The scheme works as follows: the probe and light shift laser beams propagate together through the atomic sample. 
The light shift laser shifts energetically the atomic hyperfine ground sublevels by an amount proportional to the 
intensity of the laser and inversely with the detuning of the laser from resonance. Because the light shift laser is 
tuned near the atomic Di line, it mainly shifts the ground hyperfine levels, and hyperfine components of the excited 
5p 2 Pi /2 level. The excited hyperfine levels associated with the D 2 transition are only weakly influenced by the very 
far off resonance light shift laser. However, the transmission of the probe beam is strongly influence by the shift of 
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the ground hyperfine sublevel; the effect is quite nonlinear, for the change in intensity appears as a result of the laser- 
induced shift of the atomic resonance; this quantity appears in the detuning dependence of the Beer’s Law exponent 
(the optical depth), and has a strong influence on the transmitted probe intensity. 

The effect is qualitatively illustrated in Fig. 6. The upper part of the figure is generated by first taking an image 
of the probe alone with a charge coupled device camera. Then a second image is made of the probe with the FORT 
present, and the two images subtracted. The FORT atoms themselves are illuminated with the magneto optical trap 
laser beams, making essentially a flash image of the sample. The darker spot represents the probe beam transmission 
profile; it is several pixels in size (the pixels are 13 ftra x 13 /rm. The additional influence of the light shift laser 
is evident from the lower panel in Fig. 6. In this image, generated in the same fashion as the top panel, but with 
the light shift laser present, the dark spot is unresolved at the single pixel level. We estimate the spot size to be 
several microns, consistent with model estimates of the influence of the light shift laser on the probe transmission. 
Additional details inclu ding measurements of the shift and spectral profile of the Fq = 2 —> F = 3 hyperfine transition 
are presented elsewhere [30(. Further development of this approach to creating and controlling a soft cavity in a dense 
atomic cloud is underway. Possible applications include study of Anderson localization of light in a cold atomic gas 
and studies of optical gain and the influence of atomistic disorder in quasi-one-dimension. 


Concluding remarks 

In this paper we have shown an example of the conditions where a spontaneous Raman process can reach an 
instability point associated with the threshold of random laser generation. At present there is no microscopic quantum 
theory of the random laser above threshold and in the saturation regime and the performed discussion partly explain 
that such theory would be not so easy to develop. There are mostly phenomenological approaches presented in the 
literature: some extension of a diffusion model is done in [32| . and for more information and representative references 
we address the reader to reviews fl8|-[20l| . The coherent and correlation properties of such radiation is a most intriguing 
problem. 

As known from the standard theory of a conventional cavity laser j33j , the laser radiation below threshold performs 
an amplified spontaneous emission of the atoms of an active medium considered as a random Langevin source. The 
spectral properties of the initial spontaneous emission are further modified by the amplification gain and by decay 
dynamics of the cavity mode. In this sense we see a certain analogy of the scheme displayed in Fig.[4]and the standard 
single mode subthreshold laser generation in a cavity. The source term performs the spontaneous emission, which 
would exist in any type of pumping mechanism in any laser and we emphasize that just amplification of this weak 
radiation source should be associated with a precursor of further laser generation. The main difference between cavity 
trapping and soft cavity trapping is that the latter does not select any specific spatial mode and the spectral properties 
of the emerging light are mostly controllable by the gain spectral profile. 

Considering the process above threshold it is evident that the independent multiple scattering approach is not self- 
consistent to describe the feedback and wave nature of stimulated emission. The above kinetic-type approach can only 
track the amplification process and be applicable to identification of a steady state energy balance between the active 
medium and the emitted radiation. The precise theory should turn us back to the complete diagram expansion of the 
correlation function and to the entire Green’s functions formalism substantially extended up to the saturation regime. 
In addition, we should also introduce a self-consistent master equation for the density matrix of atomic subsystem or 
more strictly for the atomic correlation functions in the complete form performed by diagram equation. 

In the considered example of the Raman laser the source generation originally occurs only in a narrow spectral 
domain in the radiation spectrum associated with the rate of spontaneous Raman emission and in a spatial directions 
along a cylindrical cavity shown in Fig. 01 If the amplification process, described by a Bethe-Salpeter-type equation 
(10.21) . involves high scattering orders then it would become more effective for those modes, which are closer to 
the Raman resonance conditions. That would make the output radiation more monochromatic and therefore more 
coherent. We also expect that spatial inhomogeneity is an important requirement, that could lead to the spatial self¬ 
organization of the emission in the saturation regime. So we point out the difference between the standard radiation 
trapping phenomenon and random lasing effect as well as between conventional laser and random laser. The quantum 
nature of the state of light also seems a challenging problem to think about and potentially can turn us to novel 
sources of quantum light, i.e. an ordered single photon regime, correlated photon pairs, etc. 
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F=4 



Fig 1: (Color online) Energy levels and excitation diagram for the Raman process initiated in an ensemble of 85 Rb via the 
Fo = 2 —> F = 3, 2 —> Fo = 3 transitions. The Raman emission is a result of the simultaneous action of microwave and 
optical u> c excitations, and both are linearly polarized along the quantization direction. The emitted light is partially trapped 
on the closed Fo = 3 —»• F = 4 transition in the optically thick atomic ensemble. While propagating through the atomic sample 
this light stimulates additional Raman emission and may increase the output fluorescence emerging from the sample. 









0 = 0 ° 


8 


k 



Fig 2: (Color online) Excitation geometry: the linearly polarized control mode (shown as a pink arrow-tube) and microwave 
field (not shown) initiate the Raman emission for the process shown in Fig. [T] The observation channels are parameterized by 
the polar angle 0 = 0°, 45°, 90° from the control mode polarization direction. 
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Fig 3: (Color online) Intensity of the light scattered at the angle 9 = 45° to the polarization direction of the driving fields, 
see Fig. [2] The frequency of the optical mode is varied in the vicinity of the Fg = 2 —► F = 4 forbidden transition and A c is 
the respective detuning. The lower panel shows the contribution of the inverse anti-Stokes scattering, the middle panel is the 
contribution of elastic scattering, and the upper panel gives the sum for both the channels. The curve thickness in the plotted 
graphs is associated with the optical thickness changed from lower to higher values bo = 1, 5,10,15, 20. 



Fig 4: (Color online) The excitation geometry, energy structure and transition diagram for the spatially inhomogeneous atomic 
system. The atoms located in the cylindrical volume crossing the middle part of the cloud have the lower energy level shifted 
by the light shifting laser operating near resonance with a hyperfine component of the Di-transition. After applying a 7r-type 
microwave pulse these atoms populate the hyperfine transition Fo = 2. The control mode is tuned in resonance with these 
atoms and initiates the initial Raman-type spontaneous emission. The emitted light is trapped by such a soft cavity and at 
certain conditions can enter the random lasing regime. 
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Scattering Orders 



Scattering Orders 


Fig 5: (Color online) Intensity distribution over the scattering orders for a pointlike dipole-type light source emitting light from 
the center of an atomic sample. The light emerges from the trapping system along the channel shown in Fig. [4] and is detected 
at infinity at the angle of the channel direction. The gray dashed curve on both the panels refers to the intensity distribution 
in different scattering orders without the optical driving field, i.e., without amplification. Other curves show the Raman-type 
amplification induced by the control field (with the reduced transition matrix element V [2f|) of different pump intensities with 
the Rabi frequencies 2V = 3 O 7 (green curve) and 2V = 407 (red curve), where 7 is the natural decay rate. In this round of 
the Monte-Carlo simulations we took into consideration only the incoherent scattering mechanism and ignored in our estimate 
additional contributions to the trapping process due to coherent scattering from the channel boundaries. Nevertheless the left- 
and right-hand panels, plotted for slightly different values of the sample optical depth b 0 = 30 and bo = 50, clearly indicate 
that the amplification process approaches unstable behavior for a critical optical depth. This instability can be associated with 
transforming the emission process towards the regime of random laser generation. 



Fig 6 : (Color online) The upper panel illustrates absorption of a tightly focused probe beam on the Fo = 2 —>• F = 3 transition 
of 87 Rb, as measured in forward scattering. The lower panel shows the influence of a light shift laser on the process; the 
noticeably smaller dark spot, which shows the modified attenuation of the probe beam, is several microns in diameter. 






































